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In this article, we derive the formula for the Casimir force acting on a piston made of real material 
moving inside a perfectly conducting rectangular box. It is shown that by taking suitable limits, one 
recovers the formula for the Casimir force acting on a perfectly conducting piston or an infinitely 
permeable piston. Lifshitz formula for finite temperature Casimir force acting on parallel plates 
made of real materials is re-derived by considering the five-layer model in the context of piston 
approach. It is observed that the divergences of the Casimir force will only cancel under certain 
conditions, for example, when the regions separated by the plates are filled with media of the same 
refractive index. 
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00 ■ I. INTRODUCTION 

Casimir effect [l| has wide applications in different areas of physics such as quantum field theory, gravitation and 
. cosmology, atomic physics, condensed matter, nanotechnology and mathematical physics (see e.g. |2Hl2j). According 
I ' to Casimir, there is an attractive force 

^: 

^SJ . acting between two perfectly conducting plates with area A and separation distance a. Naively the Casimir energy 
^ ' due to a quantum field is the sum of the zero point energies of the quantum field: 

CNi : 

G^ i Eci^^-^uj, (2) 

rn ; 

QQ , where the sum runs through all eigenfrequencies uj of the field. Generically, this sum is an infinite sum, and Casimir 
is the first one to extract a finite force ([l} acting between two parallel perfectly conducting plates. Since then, various 
methods have been proposed for the regularization of the Casimir energy ([2]). 
■ Casimir force is a pure quantum effect which can be considered as the macroscopic manifestation of the London- van 
^ [ der Waals force between atoms and molecules. Lifshitz [Tsl - fTH has laid a theoretical foundation for the description 
• of both the London- van der Waals force and the Casimir force. The celebrated Lifshitz formula expresses the finite 
temperature Casimir force between two plane parallel plates as some functionals of the electric permittivities of the 
^ . plate materials. This formula has been re-derived using different approaches ^^'^ main tool for the 

■ - - ' theoretical investigation of the Casimir effect on real materials at nonzero temperature. The original Lifshitz formula 
models the two parallel plates as two semi-infinite slabs which have infinite thickness. This approach is justified if 
the skin-depth of the material of the two plates is much smaller than the thickness of the plates. A refinement of 
the Lifshitz formula where the finite thickness of the plates is taken into consideration has been derived and usually 
related to the consideration of the multi-layer models [l6l-[23' 



Recently the issue of divergences in Casimir effect was discussed extensively i n l24|-|27ll . In connection to the 



application of Casimir effect in quantum field theory, this has been pointed out in [28|, |29| . Recently, this issue was 
reconsidered in [30l - l33| where it was argued that surface divergence terms in the Casimir energy cannot be removed 
by any renormalization of the physical parameters. Although the arguments of [30l-l33| have been refuted in the works 
[3^ [Saj. one still need to be cautious in regularizing the Casimir energy. In [sd]. Cavalcanti introduced a geometric 
setup called piston and showed that the Casimir force acting on a piston is free of divergences, due to the cancelations 
of the divergences from the two regions separated by the piston. By taking certain limits, one can recover the Casimir 
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force acting on a pair of parallel plates. Therefore, the piston approach provides an interpretation for the formula of 
the finite Casimir force acting between parallel plates. Since then, Casimir effect in the piston setup has attracted 
considerable interest [37l - l45l |. It has also been used to explore the Casimir effect in the presence of extra dimensions 
(46l - l53j . Nevertheless, these works either considered scalar fields with Dirichlet or Neumann boundary conditions 
or more general Robin boundary conditions, or electromagnetic fields with perfect electric conductor conditions or 
perfect magnetic conductor conditions. In these cases, the piston becomes a barrier for the quantum fields to penetrate 
from one region to the other, and therefore one can compute the Casimir energies in the two regions separated by the 
piston independently. The divergent part of the Casimir energies in the two regions would give rise to Casimir force 
of opposite direction on the piston and therefore cancel each other. 

Due to the success of the piston approach in the parallel plate scenario, one would expect a similar approach can 
be used to derive the Lifshitz formula for the Casimir effect on real materials. As a matter of fact, piston type 
approach has been employed by Schram [s^ j in 1973 as a regularization scheme to derive the Lifshitz formula at 
zero temperature. The generalization to include the thermal effect is still lacking. Moreover, in the work of Schram, 
some subtraction scheme was employed to get rid of the divergences, instead of observing an automatic cancelation of 
divergence as in the recent works on pistons. In this work, we re-derive the Lifshitz formula using the piston approach. 
In fact, we consider a more general setup. First we derive the Casimir energy for a three-layer model, i.e., three parallel 
layers of different materials with electric permittivities £1,62, £3 and magnetic permeabilities fj,i, fj,2, fJ-a- Using the 
piston approach, we first embed the layers inside a large rectangular box which has perfectly conducting walls serving 
as impenetrable barriers to the electromagnetic field. The finite temperature Casimir energy for the electromagnetic 
field inside the box is computed using mode sum approach. An exponential cut-off is used to regularize the Casimir 
energy. The result is written as a divergent part (which goes to infinity if the cut-off parameter approaches zero) plus a 
finite interaction term (which is independent of the cut-off parameter). Treating the second layer as a non-deformable 
piston, we compute the Casimir force acting on a piston made of real material by differentiating with respect to the 
distance between the piston and one end of the rectangular box. The limit where the other end of the rectangular 
box is brought to infinity gives the Casimir force acting between a perfectly conducting plate and a plate made of 
real material. It is shown that the divergent part of the Casimir force will only be canceled under certain conditions, 
for example, when the regions separated by the piston or the plates are filled with media of the same refractive index 
n = c^yeJI. In general, one cannot ignore the cut-off parameter. 

To derive the Casimir force acting on two parallel plates both made with real materials, it is necessary to consider 
the five-layer model. Using the same approach as the three-layer model, we derive the Casimir energy and the Casimir 
force for the five-layer model. In the limit where both ends of the confining rectangular box go to infinity, we obtain 
the Casimir force acting on two parallel plates, both made with real materials. The limit where the thickness of the 
two plates are infinite recovers the celebrated Lifshitz formula. As in the case of the three-layer model, we observe 
that the divergences of the Casimir force would cancel each other if the three regions separated by the two plates are 
filled with media of the same refractive index. 

Finally, we would like to remark that in this work, the piston approach does not really refer to the piston, but 
rather the finite rectangular box with perfectly conducting walls that serve to confine the field into a finite region. 
The advantage of this approach is that in solving for the field modes, one do not have to use imaginary energy modes 
as in some of the previous works. This bridge the gap between the derivation of the Casimir force for real materials 
and the derivation of the Casimir force for ideal metals. 



We consider a piston with thickness d2 = b — a dividing a rectangular box of dimension L x L2 x ^3 (see FIG.[T]) into 
two regions filled with materials with electric permittivities and magnetic permeabilities £1, fii and £3, ^3 respectively. 
The piston is assumed to have electric permittivity £2 and magnetic permeability /i2 and the walls of the rectangular 
box is assumed to be perfectly conducting. Maxwell's equations take the usual form 



II. CASIMIR ENERGY INSIDE A PISTON SYSTEM 



V.D ^pf, 



V.B =0, 




(3) 



where D = £E and B = /_tH. As usual, we introduce the potential $ and A such that 




B = V X A, 
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FIG. 1: A piston with electric permittivity £2 and magnetic permeability fj,2 diving the rectangular box into two regions, 
which are filled with materials having electric permittivities £1,63 and magnetic permeabilities /il,/^3. 



and impose the radiation gauge 



V.A = 0. 



In the absence of free charges and free current, i.e., Pf = and J f — 0, we can impose the additional condition $ = 0. 
Then the first three equations in Q are automatically satisfied and the fourth equation give 



5' v.^ 



A = 0. 



(4) 



For the boundary conditions, the continuity of D . n, E x n, B . n and H x n across the boundaries implies that eEi, 
E2, E3 and Bi, —B2, --B3 must be continuous in the xi direction. On the other hand, the perfectly conducting 
conditions on the boundary of the rectangular box imply that on the boundaries Xi = and xi — L, E2 = — 
and Bi ~ 0. On the boundaries X2 = and X2 — L2, Ei = E^ = and B2 — 0. On the boundaries X3 = and 
X3 = L3, Ei= E2^0 and B3 = 0. 

A solution of equation ^ satisfying the gauge and the boundary conditions can be written in the form 
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'k2' 
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. 7rfc2X2 . TTk^Xs 

(xi) sm — sm — e 



L 



k2 d(f){xi) TTk2X2 . nk^xs 
A2 —- cos — - — sm — e 
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L2 dxi 
kz d4){xi 



L2 
nk2X2 



Az —- X sm — cos — e 



L3 



dxi 



(fc2,fc3) e^TM :-N2, 



for TM modes, and 

Ai =0, 



. h If . Trk2X2 . TTksXs 

M =— cos — sm — - — e , 

L-j, L2 L3 



A. 



k2 , nk2X2 irksxa 

— V-'(2;i)sm cos e , 

L2 L2 L3 



(fc2,/c3) e^TE :Ng\{(0,0)}, 



for TE modes. Here No = N U {0}. The functions £(^(a;i) and 



d<t>{xi) 



are continuous, and (t>{xi) satisfies 



= 0, 
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where p{xi) = pj for xi in region j whose electric permittivity and magnetic permeability are Sj{uj) and fJ-ji^^) 
respectively. The functions ip(xi), ~ ^^q!^^^ continuous, and ^(xi) satisfies 

d^i^ixi) 



Finally, for j = 1,2,3, 



dx\ 



ej{ijj)^j{uj)oJ^ 



Xl), 




^(0) 


7rfc2 


2 




[ L2 \ 


+ 


^3 



For the fimction (j){xi), assume that 

4>{xi) 



Age^P^"^! + Age-^P^^i, 



< < a, 
a £ a^i < 6, 
6 < xi < L. 



Then the boundary conditions give six equations for the six variables Ai, . . . , Ag: 

Ai - A2 = 0, 

ei (Aie'Pi" + Aae-^Pi") = £2 (Age^P^" + Ue-'P^") , 
ipi (Aie'Pi" - K2e-"P''') = ip2 (Aae'^^" - A^e-'P'") , 
£2 (Age^P^^ + A^e-'P^'') = £3 (Age'P^'^ + Age-^f^'') , 

, se^P^'-Aee- 

56"— - Age-"^^ " 



(5) 



These homogeneous system in (Ai, . . . , Ag)^ has a nontrivial solution if and only if the determinant of the corre- 
sponding matrix is zero. However, we can also solve the system in the following way. The first equation and the last 
equation give 



A2=Ai, As^Ag 
The second and third, fourth and fifth equations give respectively 



1 
2 
1 

2 \-r. 



~i{pi+P2)a 



pi{pi+P2)a „+ „-iipi-p2)a 



12 
12 

^ p-i(P3+P2)b 



32 



Ai 
A2 
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Ag 



+ pi{Pl-P2)o- 



''126 



^12^ 



i(pi+P2)a 



-i{pi+P2)a 



+ ^-i{pi-P2)a ] Al) 



' 12"^ " ■ - -r ^126 
T.+ pi(P3-P2)b „-2ip3L 

2 V-r7ne*(P3+P2)''e"^^P3-^ +rt,e~*(P3-P2)6 ) «' 



-^p-i(p3+P2)6 



' 32^ 



(6) 



(7) 



where 

^ Pk £fc 

Equating the two equations in ([7]) give an equation of the form 



Al 



Ag 



which has a nontrivial solution if and only if aS — = 0. In other words, the system ^ has a nontrivial solution if 
and only if 



- p-ipidi^ ^^+^ip3d3 _ rr_p-iP3d3\ ^ip2d: 



^12^ 



^32 6 



0, 
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where di = a, d2 = b — a and = L — b. Discarding the nonzero factor e ^p^^ and define 

- (r+ e-'fi'^i - r^^e'P"'') (r+ e-^f^"^^ - ^ale'"^'^^) e-'P^''^ fc = (^2,^3) € 5tm 

where for j = 1, 2, 3, 



The set of TM energy eigenmodes is the union of the sets of real nonnegative zeros of Fo;tm('^, k), k = (^2, ^3) € 5tm- 
For the TE modes, one can show analogously that the set of TE energy eigenmodes is the union of the sets of real 
nonnegative zeros of 



(9) 



where 



Compare © to (|H]), we find that Jo;TE(i^,fe) is obtained from ^o;TM(w,fc) by replacing by s^. In the following, 
we assume that all the zeros of Fo;tm('^, k), k — (^2, fc3) € 5tm, and /o;Te('^, k), k = (^2, k^) e 5te, are real. 
The finite temperature Casimir energy inside the piston is given by 

EcM - ^ E + ''bT l°g (1 - (-^)) ' (10) 

modes modes 

where we have introduced a cut-off parameter tc to render the first sum finite. Using the generalized Abel-Plana 
summation formula (ssl. |56| (see Appendix IX)) . we compute this finite temperature Casimir energy in Appendix IB] 
Up to the terms that give non-trivial limits when tc ^ 0^ , the result is 

3 

£^Cas(tc) = dj^ji^c) + A£^Cas, (11) 



where 



5,(tc) = -f E + E 1 / |^e-- + fcBriogfl-expf--^))jdp,(..,fe), (12) 



and 

knT 



Ai?cas =^ E E (l - A™(z6,fc)e-^''^(«'^'=)'^^) (l - A™(z6,fc)e-^''^(«'''^)^^; 



(13) 



with 



2TT\l\kBT 
I = ^ : 



^■'^^'^ e,(*Og,(e,fc)+£,(zOft(e,fc)' ^ M,(*e)%(e,fe) + Mfe(*e)<z,(e,fe)' 

Notice that X]j=i djE,j{tc) contains all the tc — > 0+ divergences. The term di'E.i{tc) is independent of d2, d^ and the 
parameters of media 2 and 3. It represents the self energy of media 1. The term Si(tc) can be interpreted as the Casimir 
energy per unit length that would exists in the region between xi = and xi = a if the boundaries at zi = and 



6 



xi = a were absent. Similar interpretation can be given to the terms d2S2(ic) and d^'E.^ltc). The term AEcas is called 
the interaction term. In the literatures on Casimir effect, AE'cas is usually regarded as the regularized Casimir energy 
by claiming that the term dj'B.jitc) should be subtracted away as the energy in the absence of boundary. This 

regularization scheme is a little superficial, since we have to consider the three regions that have different properties 
separately. As explained in [2J-|27[, one should not ignore the divergent terms, but should instead give a physical 
interpretation to the cut-off parameter tc- Nevertheless, in the special case Si(ic) = '^2{tc) — S3(tc), then the term 
diE.i{tc)+d2^2{tc)-'rd^'B.^{tc) = LSi(ic) is independent of di, ^2, c^s, and the subtraction of this term from the Casimir 
energy becomes natural since it would not contribute to the Casimir force (see next section). Si(ic) = '^2{tc) = '^sitc) 
holds when for example, the system is isorefractive, i.e., the refractive index n(uj) ~ c^J e{uj)^{uj) is the same for the 
three materials. 

The interaction term of the energy AE'cas shows that in the high temperature limit, the leading term of the energy 
is linear in T, given by the sum of the / = terms. The zero temperature Casimir energy is obtained by taking the 
limit T -> in ^ and dTS]), which gives 



where the interaction term is 



4 ^ /°°log|(l-A™(*f,fe)e-^^^«-'=)'^^)(l-A™(*f,fc)e-2^3(«..).3j 

- (e-2«^(«''=)'^^ - A™(zC,fc)) (e-2«3(«''^''^3 - A™(zf,fc)) e-^'J^K.fe)^^ + (tM TE) . 

In the limit where the two ends of the rectangular box a;i = and xi — L are brought to infinity, i.e., di oo and 
da — > oo, the interaction term of the Casimir energy (I13p gives 



Ai?cas=^ E E l°g(l-^™(*^''^)^™(*^''^)^"''^*''''^'^)+(TM-^TE) 



In addition, if the media are infinite in the transversal direction, i.e., L2, — !■ 00, then A| has to be replaced by fc^ 
and the summation over k has to be changed to an appropriate integral. More precisely, in the limit L2,L3 — >■ 00, 
the interaction term of the Casimir energy density is 



Afcas = ^f:f {log (1 - A™(*6,fc)A™(z6,fc)e-2^^(«-'=)'^=) + (TM TE)} kdk, (15) 

where qj{£^,k) — \/sj{i(,)fij{i(,)(,'^ + k"^. This is precisely the Lifshitz formula [13] for the Casimir energy density 
when two semi-infinite media with permittivities and permeabilities £1,^1 and £3, /is are separated by a medium of 
permittivity £2 and permeability /i2. Notice that we have ignored (or subtracted away) the divergent part of the 
Casimir energy c?iSi(ic) + d2'^2{tc) + d^'^j,{tc) and only take the limit di, ^3 — ?> cxd on the interaction term. Later on, 
we are going to show some situations where the divergences would naturally disappear when we consider the Casimir 
force. 



III. CASIMIR FORCE ACTING ON THE PISTON 



We assume that the piston is allowed to move freely inside the rectangular box, but it is non-deformable, i.e., its 
thickness d2 = h — a is fixed. The position of the piston can be described by the variable a = di. As a function of a, 
di = a, d2 is fixed and ^3 = L — a — d2- The Casimir force acting on the piston is given by 

ppiston/ X _ dEcas _ dEcas dEpas 

da " ddi dds ■ 

From (HI]), dm) and ([13]), we find that 

F^ir (a) = S3(i.) - + AF^::°^{a), 
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FIG. 2: The Casimir force Jcas*°"('^) ^ function of the plate separation a for different values of £2,r ~ £2/£o and 
M2,T- = /i2//io- Here L = 2cm, L2 = I/3 = 1cm, d2 = 1mm and T — IK. The two regions separated by the piston are vacuum. 
For the graph on the left, /i2 = Ho- For the graph on the right, £2 = £0. 



where Si(ic) and 'E.^itc) are defined by (|12p . and the interaction term is 



with 



Ls-TM^-ksT ^ ^ (1 - A™e-2gidi) (1 _ A™e-293d3) _ (atm „ g-29idi) (a™ - 6-293^3) £-292^2 ' 



A-Fcas-TM obtained from '-^F^i^.rj.j^ by interchanging the index 1 with 3, and ^F^l^r^.^ is obtained from ^F^l^r^,^ 
by changing A™ to A™. Unlike the piston cases that have been considered so far, there are still divergences 
^3{tc) — Si(tc) in the Casimir force coming from the self energy terms. This term will vanish if and only if !E!i(<c) = 
'^sitc), which would happen for example when the two media separated by the piston is isorefr active. In particular, 
if the two media 1 and 3 are both vacuum, which is the case in the previous works that have considered the piston 
scenarios, then there are no divergences, and regularization is not required for the Casimir force. 

For the particular case where ei = £3, /i = /ia, the self-energy terms Si(tc) and '^.^{tc) cancel and we can set the 
cut-off parameter to zero. Let q — qi — qs and A"'"'^/™ = A™^"""^ = A™'''^^. Then the Casimir force acting on 
the piston is given by 

ppisto„^_, 7. V V gA^M (,-2,.._g-2,.3)(l„e-2^2d2) 

Cas n _ A™e-2<Z<il) (1 _ ATMg-2gd3) _ (ATM _ f,-2qdi\ l/\TM „ g-2gd3) g-2g2<i2 

^GiStm I— — 00 V / 

+ (TM — > TE) . 

From this, we find that the sign of the Casimir force is determined by the sign of A"^'^, A™, and the sign of di — d-^. 
Notice that in general the Casimir force depends on the thickness of the piston, and it vanishes as the thickness of 
the piston goes to zero, unless A™ = ±1 and A"'^^ = ±1. FIG. [Hand FIG.[3]show the special case where the regions 
separated by the piston are vacuum, i.e., ei = £3 = £0, /^i = ("3 = A*o and the piston is made of material with constant 
permittivity £2 and permeability //2- We find that when £2 > £0, the Casimir force attracts the piston to the closer 
wall, and when /i2 > /xg, the Casimir force push the piston away towards the middle. This can also be directly verified 
from (IT7)) . 

In the special case where the regions separated by the piston are vacuum, i.e., £1 = £3 = £0, /^i = /i3 = /io, and the 
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FIG. 3: The Casimir force F^^°^(a) as a function of the piston thickness 6.2 for different values of S2,r = E2/E0 and 
A'2,r = A'2/A'o. Here L = 5cm, L2 = L3 = 1cm, di = 5mm, T — IK. The two regions separated by the piston are vacuum. For 

the graph on the left, ^2 = Mo- For the graph on the right, £2 = eo- 



piston is perfectly conducting, achieved by the limit 62 ^ 00, we find that 



ATM... .X ATM... ^2\J^ + Xl- Soy/e2H2^'f + 



£2 V SoV£2f^2^'t+^i 



(18) 



When I = 0, the limit of the second term has to be taken using the prescription of Schwinger, DeRaad and Milton 
[stI where one first takes the limit £2 — ^ 00 before setting I = 0. Alternatively, one can also take the limit 112 — >■ 0. 
In fact, superconductors can be considered as perfectly conducting and it is well-known that for superconductors, the 
magnetic permeability is zero. Substituting the limits ([T5)) into the Casimir force (|17p . we find that 



(19) 

The two terms in the brackets represent the Casimir energies from the left region (with thickness di = a) and from 
the right region (with thickness d2 = L — a — d2) respectively. This formula is exactly the same as the formula for a 
perfectly conducting piston moving freely inside a perfectly conducting rectangular cavity derived in [sol l43j . Notice 
that this formula is independent of the thickness of the piston d2. 

In the case the piston is infinitely permeable, the limit is obtained by /12 — ^ 00. As dual to the perfectly conducting 
case, we stipulate that £2 ^ 0. Then 

-1, A™(zC,, fc) -1. 

In this limit, the Casimir force (|17p becomes 



^piston ^ j^^rj. ( ^ + ^ ) ^ ^ V V he ) ~^ V I he ) ^ 



kes^^ kes^J i=-^\exp(2dJ{'-^Y + Xl)+1 exp 2^ V (^^)' + A| + 1 



(20) 
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which coincides with the result for an infinitely permeable piston moving freely inside a perfectly conducting piston 
we derived in |4J1 . Notice again that the Casimir force is independent of the thickness of the piston c?2 • 

As one can see, in general the nature and the strength of the interaction term of the Casimir force acting on the 
piston depend on the properties of the materials and the thickness of each material. It does not decouple into 
the Casimir force acting from the left region and the Casimir force acting from the right region. In the limit the right 
end of the rectangular box is moved to infinity, i.e., da — >■ oo, we find that A_F'(^jjg.rpj^ — >• 0, Ai^(?jjg.rpg — > 0, and 

AF«-- = -ksT E E '^M \atm'\tm,^., J -A + (TM ^ TE) . (21) 

This formula can be interpreted as the interaction term of the Casimir force acting on a perfectly conducting plate 
and a plate made of real material (with permittivity £2 and permeability ^2) which are embedded orthogonally in an 
infinitely long rectangular cylinder, with the medium between them having permittivity ei and permeability fix and 
the medium outside them having permittivity £3 and permeability (13. (j2ip shows that when the separation between 
the plates di gets large, the magnitude of the interaction term of the Casimir force decreases exponentially. 



IV. CASIMIR FORCE BETWEEN TWO PERFECTLY CONDUCTING PLATES SEPARATED BY 

MEDIUM WITH CONSTANT REFRACTIVE INDEX 

If both plates are perfectly conducting so that A™ — A™ — A™ — A™ ~ 1 , and the medium 1 between the 
plates has constant refractive index ni — c-y/£i/ii, then we find that the interaction term of the Casimir force (j2ip 
becomes 



AFc,s= 2- E E 7 / 2 =\ ■ (22) 

(fe2,fe3)ez2\{o}i=-<x) exp ( 2diy^n2 (2^^^) + AM - 1 

This term is always attractive, and its magnitude gets smaller if the medium between the perfectly conducting plates 
has higher refractive index. In the zero temperature limit, the interaction term of the Casimir force (j22p is 



00 




E 

ki=i (fc2,fe3)ez^\{o} 

which is ni times smaller than the corresponding Casimir force on perfectly conducting plates separated by vacuum. 
When the plate separation di is much smaller than the size of the cross section, i.e., di <C -^2,^3, the leading terms 
of the interaction term of the Casimir force is 



hcTT^L2L3 

240nidf ' 



Afcas=- ,4 ' ifT = 0, 



4^, ' ^id,T»l. 

Notice that when T — 0, the leading term is ni times smaller than the corresponding term for perfectly conducting 
plates separated by vacuum, but when the temperature is large enough, the leading term becomes independent of the 
refractive index of the medium between the plates. 

The interaction term has not taken into full account the influence of the material between the plates. There 
is an additional contribution to the Casimir force arising from the difference of materials in the region between the 
plates and the region outside the plates, given by 

S3(ie)-Si(ie)--( E + E 1 / (^e-*^" + fcsTlogfl-expf--^)Hdp3(a^,fc)-(P3^Pi). 



10 




FIG. 4: A five-layer model 



Assume that the permittivity 63 and the permeabihty fi^ of the medium in region 3 are also constants with refractive 
index n^. Then up to the term constant in tc, 



Sl(ie) 



3fe^2^3 3,-4 
- ^2^3 "l^c 



47rc 



' = 1 (fc2,fe3)eZ2\{0} 



Cfl(3) 



+ 



1 



00 00 



2. 



720 L3 167rL2 ' ^ ^ 

^ -i fc2 = lfe3 = l 



fe) +(&)• 



/ 27rfc2fc3L3 




(23) 



and S3(ic) is obtained from this formula by replacing ni with 77.3. It is obvious that if ni 7^ 713, the terms Si(ic) and 
^sitc) due to the self energy of the system are not equal. As explained in [23 - [27j . in this case, the parameter tc has 
to be present in the expressions for Casimir energy and Casimir force. 



THE FIVE-LAYER MODEL 



The piston system we consider in Sections and IIIII can be considered as a regularization device for a three-layer 
model. We have asserted that in the limit where one end of the rectangular box is brought to infinity, we obtain the 
Casimir force acting on a pair of parallel plates, one being perfectly conducting and the other made of real material. 
In this section, we consider the more general case of the five-layer model (FIG.S]) which is the more appropriate model 
for two parallel plates both made of real materials. 

The Casimir energy and Casimir force of a five-layer model have been derived by some authors using different 
methods [TB-Hl] • We are going to use the piston approach. Namely, first we assume that the whole system is inside 
a large perfectly conducting rectangular cavity [0, L] x [0, L2] x [0, L^]. As depicted in FIG. SI there are five layers of 
materials with electric permittivities ei, £2, £3, £4, £5 and magnetic permeabilities /ii, /i2, /i3, /is from left to right. 
The first medium extends from xi = to =01, the second from xi — ai to xi =02, the third from xi ~ 02 to 
xi — a^, the fourth from xi = 03 to xi ~ 04 and the last one from xi — 04 to xi — L. 

Let di — fli, (i2 = 02 — oi, c?3 = 03 ~ 02, (i4 = 04 — 03 and = L — 04. Then using the same method as Section HIl 
we find that the Casimir energy of the system is 



Casj 



(24) 
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where ^j{tc) is as in (jl2[) . and the interaction term of the Casimir energy is 



[1 + A™A™e-2'?4d4] _ [^TM ^ Ajg^e^^g.d,] ^-2,5^5^ _ ([g-^^^"^^ + A™ A™] e-^^^'^^ - [A™e-2?2d2 ^tm] j 



X 



([e"2g4d4 + A™ A™] e-295rf5 _ [ATfe-^?^"^^ + A™]^-^'''^' I + (TM ^ TE) 



where q, = q,{^i,k), A™ = A™(*6, fe), A™ = A™(*6, fe). 

The functions qj{^i,k), A™ = A™(i6,fc) and A™ = Aj^ {i^i,k) are defined as in (0. 

Now assume that medium 2 and medium 4 are two non-deformable pistons/plates embedded orthogonally inside a 
rectangular box, dividing it into three regions filled with different materials. The Casimir forces acting on the plate 
on the left (medium 2) and the plate on the right (medium 4) are given respectively by 

- " dd, + dd, - -^'^^'^ + ^'^^'> ~ ^dT ^ ~ddr' 
_ dEcM . dEcM „ . „ , dAEc^s , dAEc.^ 

dd, + dd, - + ^'^'^^ - -ddT ^ ~^dr- 

By moving the two auxiliary ends of the rectangular box to infinity, i.e., di,dQ — )■ 00, since 

dAEcus rfi^oo^ Q dAEcas d^-^oo^ p 

ddi ' 9^5 

we find that 

F^,, --Si(t,)+S3(U+AFcas, 
J^Cas=-23(y+S5(ic)-AFcas, 

where the interaction term of the Casimir force acting on the two plates is 
dAEc^, 



AFcrx — ., , , _ 




(25) 



A™) (ATMe-294d4 + A™) 



The divergences — Si(tc) + S3(tc) and —'^^{tc) + '^sit^) come from the self energy of the system and in general do 
not vanish. A special case they would vanish is that the media 1, 3 and 5 are isorefractive. In the symmetric setup, 
i.e., media 1 and 5 are isorefractive, then the force —Ei{tc) + '^sitc) acting on the plate on the left and the force 
— S3(tc) + 55(^2) acting on the plate on the right have opposite signs. In the limit the plate on the left is perfectly 
conducting, i.e., £2 —J' 00, then A™, A™ — > 1 and A™, A™ — >■ —1, we find that the formula (|25p coincides with 
the formula (|21l) after the re-indexing 3 ^ 1,4 — 2 and 5 — ?► 3. This is consistent with our interpretation of formula 
(j2ip as the Casimir force acting between a perfectly conducting plate and a plate made of real material. 

Usually in the works of Casimir effect on parallel plates, it is assumed that the size of the cross section of the plates 
is much larger than the separation between the plates. In this limit, the interaction term of the Casimir force density 
acting on the plates is 



kgT ^ / l + A2^™A|™e-^«2''2 l + A!™A,Y'e-^''"'M , . \ 

(26) 

In the zero temperature limit, (j25p gives 



-1 



AF--- f °" , f (1 + A?3^A7,^e-^^--'-) (1 + A,TAi-M,-2,...) de + (TM TE) 



(27) 
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FIG. 5: The Casimir force -Fcas as a function of the plate separation da (left) and the plate thickness d2 = (right). Here 

L2 — Ls = 1cm, T — O.IK, the regions separated by the plates are vacuum. For the plate on the left £2 ~ 2eo and 
/i2 ~ /i2,r/io, where fi2,r = 3.6, 3.8, 4.0. For the plate on the right, £4 = 4eo, A** = f^o- For the graph on the left, the thickness 
of the plates are d2 — d4 = 0.5mm. For the graph on the right, the separation between the plates is da = 1mm. 



In addition, if <C L2, L^, then the zero temperature Casimir force density acting on the plates is 




1 I ATM ATM -2^^2£^2^ (-1 I ATM ATM -2g4d4\ \ ^ 1 



(28) 

Formulas (|25|) . (|26|) . (f27l) and (f28|) can be used to study the nature and properties of the Casimir force acting on 
two plates made of real materials. We would like to emphasize that in the case the materials in the regions 1, 3, 5 
separated by the plates are not isorefr active, one also has to take into account the forces coming from the terms Si, S3 
and S5 and the cut-off parameter tc cannot be set to zero. 

In the case the three media 1, 3, 5 are the same, i.e., £1 = £3 = £5, and fii = /i3 = /i5, we find that the self energy 
terms cancel and the Casimir force acting on the plates is ((25|): 

In addition, if the two plates are made of the same material, i.e., £2 — £4 and /i2 = /i4, (|29p becomes 

Fcas = - fcBT ^ -^"''''^ -A + (TM ^ TE) . (30) 

.e^M^= c» V [Aj3^]'(l-e-2«..^)(l-e-2,.rf0 ) 

It is easy to see that in this case, the Casimir force is always attractive. This is a special case of the theorem [s^ [s^ 
which asserts that the Casimir force between two bodies with the same property is attractive. In general shows 
that the dependence of the Casimir force on the thickness of the plates and the distance between the plates can be 
quite complicated. In FIG. [SJ we plot the dependence of the Casimir force on the distance between the plates and 
the thickness of the plates. The graphs show that with £2 = 2£o, £4 = 4£o, ^4 — fJ-o and ^2 — 3.6/io, 3.8/xo, 4.0/zo, the 
Casimir force can be attractive or repulsive depending on the distance between the plates and thickness of the plates. 

Finally we see how we can recover the Lifshitz formula from (1261) . As mentioned before, when the regions separated 
by the two plates are filled with isorefractive media, then the divergences of the Casimir force cancel. If in addition. 
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we assume that both plates have infinite thickness, i.e., ^2,^4 00, ([26)1 gives 

Fcas = - ^ E r ([^™(*^'' fc)] [A4T(*6, fc)] e2''3(6.fc)<^3 _ 1) kdk + (TM TE) . 

/— — 00 ^ 

This is precisely the Lifshitz formula for the Casimir force acting on two parallel plates made of real materials. Thus 
we see that the piston approach can be used to derive the Lifshitz formula. An advantage of this approach is that the 
divergence of the Casimir force is not ignored by ad hoc subtraction, but is seen to vanish due to pairwise cancelation. 

VI. CONCLUSION 

In this article, we derive the Casimir energy and Casimir force acting on a piston made of real material that moves 
freely inside a perfectly conducting rectangular box. We also use the piston approach to derive the Casimir energy 
and Casimir force acting on two parallel plates made of real materials. We use an exponential cut-off to regularize 
the Casimir energy. When the regions separated by the piston or the plates are filled with isorefractive media, we 
show that the divergences of the Casimir force cancel each other and therefore the cut-off parameter can be set to 
zero. In general, the Casimir energy contains divergent terms that depend on the cut-off parameter that cannot be 
ignored. Although we have restricted our consideration to three-layer model (for the piston) and five-layer model (for 
the parallel plates), our method can be easily generalized to Casimir effect on multi-layer models. Moreover, it is easy 
to see that there is nothing special in choosing the rectangular cavity as a space cut-off for the electromagnetic field. 
We can also use a cylinder with arbitrary cross section as a substitute. For doing so, one just has to replace the 
by an appropriate spectrum. 

For simplicity, we have illustrated the behavior of the Casimir force using non-realistic examples where the electric 
permittivities and magnetic permeabilities of the materials are assumed to be constants. In reality, the electric 
permittivities and magnetic permeabilities depend on temperature and also the frequency of the field. There are some 
models that have been used for the electric permittivities such as the plasma model and the Drude model. An issue 
that we do not discuss in this article is the leading term of the thermal correction to the Casimir force at low and high 
temperatures. This has been an active topic of research and the results often depend on the particular model chosen 
for the electric permittivities and magnetic permeabilities. In the case where the self energy terms cancel, we have 
shown that the high temperature leading term of the Casimir force is linear in temperature. In case the self energy 
terms do not cancel, further regularization may be needed to cancel the high temperature leading terms in the self 
energy that are higher than linear orders [goI]. This issue will be discussed in more detail elsewhere. 
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Appendix A: Generalized Abel-Plana Summation Formula 

The generalized Abel-Plana summation formula states that if fo{z),fi{z) and f2{z) are meromorphic functions, 
and 

\im [^\fo{x + iY)~f,{x + iY)}dx = 0, for i = 1,2, (Al) 
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FLd z = a 



Re2i = P 



FIG. 6: The integration contour for the integrals in (|A2 



then by residue theorem, 



wo{z)RcSzfo{z) - ^ wi{z)ReSzfi{z) - ^ W2{z)ReSzf2{^- 

a<Rc z<P 



a<Re z<P 
Im z>0 



a<Re z</3 
Im z<0 



^ I {Mz)-hiz))dz + ^ 

/•OO 



(/oW-/2W)ciz 

J [Mx + iy)- fi{x + iy)j ^_^dy+— J \^Mx - ly) - f2{x - ly) j 
■J— Hm / \ fi{x + ie) - f2{x - ie)\dx. 



dy 



(A2) 



Here 



wo{z) 



1, ifzeSo, 
1/2, ifzeSSo, 



1, ifzeSi, 1, iizeD2, 

wi{z)={l/2, if i«2(2) = < 1/2, if z e \ {a, /?}, 

1/4, if 2; = a or /3, I 1/4, ifz = aor^. 



where Dq = {z : a <Rez < fi} , T)i = {z : a < Re z < f3, Imz > 0} , 'D2 ^ {z : a < Rez < /3, Imz < Oj. The 
integration contours L^, , are shown in FIG. [G] This formula can be proved in the same way as in [55l[56|. 



Appendix B: Computation of the Casimir energy 

In this section, we compute the Casimir energy (fTO| inside the piston system using the generalized Abel-Plana 
summation formula (jA2p . First we consider the zero temperature part 



cue 

modes TM modes TE mode 



Recall that the TM modes are the positive zeros of the functions Fo:tm(w, k),k — (^2, k^) e 5tm ©■ Therefore we 
take a = and let /3 — )• 00 in (jA2|) . Let 



+ T.+ pWldi+ip3d3+iP2d2 
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and define 

/,.TM(w,fc) = ^e-*='^^logF,.TM(w,fc), J =0,1,2. 

One can then check that the condition (IA1[) is satisfied. With the further assumption that r^2{'^T^^) ^-nd r^2{^Tf^) 
do not have zeros, we find that /i;tm(w, k) and /2;Tm('^j k) do not have poles. Therefore the generalized Abel-Plana 
summation formula (|A2p gives 



- E / (/o(*^,fc)-/i(ze,fc))de-- Yl / (/oHe,fc)-/2He,fc))rfe 

TM modes ^ fcG^TM ° ^ fceSxM ° 

- — Y] lim / {fi{uj + ie,k) - f2{uj - ie,k))duj 

fcGSxM " 

For the right hand side, a straightforward computation gives 

fo{^t k) - fM. k) =f log/ (e^^P^i^i'^)^^ - U^P^i^W^'^ _ .^^v.^m 



V r+(zC,fc) 



'"32(*C,fe) 



where 



Similarly, one finds that 

M-^^,k) - U-rtk) =Me^*e«^iog j L - :^p^-2^v.i-m\ U - ^32He, fc) -2.,3(-.C,^).. 



-2ipi 



i(-ie,fe)£;i _ ^12( ( -2ip3(-i?,fe)d3 _ ^32 ( fc) A -2ip2 (-i^ ,fc)d: 



Under the assumption that for all ^ e M, j = 1,2,3, i^) = Sj(i^), lJ.j{~i£,) = fJ,j{i£,), we find that 

fc) - v'^,HOm.HOHO'-a| = -*y^e,(*0M,(*0C2 + A| = fc) 

and 



Therefore, 



foi-^t k) - hH^, k) =y e-*=« I logj (l - A™(^C, fe)e-2,,K,fc)., j _ ^tm^^^^ fc)e-2.3K,fc)<i. 
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For the last term in ([B]), under the assumption that £j{u!)^j{u!) is real for > 0, we find that if 
then 



lim pj {u! + ie,k) —iy^i. ~ lim Pj — le, fe) — ~ ej(uj)iij{uj)Lij'^ (Bl) 



whereas if 
then 



lim Pj {uj + ie, k) = lim pj (w — ie, k) — J ej[uj)^j{uj)ijj'^ — A| > 0. 

1S-S-0+ £->-0+ ' 



Therefore 
1 



lim / {fi(uj + ie,k) — f2(uj — ie,k))duj 
2ni e^o+ Jq 



— — lim 

An 



TrY.'^n ioe'''^dp,{uj,k) 



where the part with ej{uj)fij{uj)uj^ < has been canceled due to (|B1[) . Grouping everything together, we have 



h 

2 



[modes ^^;ee5TMj = l -'p.(".fe)>0 

- E rCcos(U):^log|(l-A™(*f,fe)e-2''^(«''=)^^)(l-A™(*f,fc)e-2''3«^^^ (B2) 
2«i«,fe)di _ A™{i^,k)) (e-293(«.fc)d3 _ A™(if,fe)) e-29=«''=)'^4^f' 



Notice that the second term in (|B2p is finite when tc — > 0+. Therefore all the tc — > 0+ divergent terms come from 
the first term in (jB2l) . which is a linear homogeneous function of di, (i2, ds- Setting = in the second term, we find 
that when tr O'^, 



2 ^ 

TM modes 



^E^.^.;TM(tc)-f E r ^pog\{l-AJ^{^^,k)e-W)d.^ [l-AJ^{^^,k)e-''^^m) (b3) 
j=i ^fceSTM-^o ^ I 



where 



Xj:TM{tc) ^ TT E / uje-*'''^dpj{uj,k 



A comparison of Fo;TE(w,fc) ® with Fo;TM(w,fc) ® shows that the expressions for the zero temperature Casimir 
energy due to the TE modes can be obtained from the expressions for the energy due to the TM modes (jB3p by 
replacing Aj^ by 



jk 



TE _ ^j7c(^^^fc) 
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Next we consider the thermal correction to the Casimir energy 

Two 



SrEcas ^ ksT ^ log I 1 - exp 



knT 



modes 

Define the zeta functions 

2N 



modes modes — oo 

It is "well-known that 



6tEc^ = -^2'(0), 



where 



Using the generalized Abel-Plana summation formula again, we find that 

3 



Co;Tm(s) =- E E^^J / ^ '^'dpj{u:,k) 

+ - Yl ^sin(..)r^«|log^(l-A™(^^,fc)e-2^^«■'=)'^^)(l-A™(^e,fc^^^^^ (B4) 



and 



Ct;tm(.)4 E EE^J P + '^^^ 

+ i E E / BinM(^^-ef)-^|log (l-A™(ze,fe)e-^^^«.'=)'^^)(l-A™(ze,fe)e-^ 
_ (^g-2,,K,fe)d, _ A™(z^,fc)) (e-2«^(«^'='''^ - A™(*^,fe)) e-2,2(«,fe)d.|rf^^ 

(B5) 

where 

. 27r|?|fcBr 

Therefore, we find from (lB4l) and (iBSl) that 



■3(«,fc)d3 



TM modes ^ \ B / / .^^ keS-rm "'0 S 

■ log|(l-A™(zei,fc)e-2'i«-'=)'^i) (l-A™(zei,fc)e-2'^«-'=)'^^') - (e-^?!^^-'^)^^ - A™(zei, fc) 



2 

1= 
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where 



Together with the zero temperature Casimir energy (IBSp . we find that the contribution to the finite temperature 
Casimir energy from the TM modes is 



Yj-TM = V / log ( 1 - exp (--^]] dpj{u},k). 



Ec.s;TM{tc)^Y.do^r.TM{tc) + ^ E E log (l-A™(zC;,fe)e-2'^«-'=)'^^)(l-A™(z6,fcK 



-2q3{ii,k)d3 



(B6) 



where ,^j;TM(ic) = ^j:TM(ic) + yj:TM- The contribution to the Casimir energy from the TE modes is obtained from 
this formula by replacing A™ with A™. 

In deriving the formula (|B6I) . we have made a few assumptions on the nature of the zeros of the functions i^o;TM('^, k), 
fo;TE(i^, fe), and the properties of the permittivity and permeability functions £j(y>)-, However, in general some 

of these assumptions need not hold. For example, there may be zeros of i^o;TM(w,fc) and -Fo;Te(w,A;) that do not 
give rise to energy eigenmodes. In such cases, we will obtain extra terms for the Casimir energy. Another possible 
failed assumption is the assumption that ej{—i£) = ej{i£,) and = this case, the second sum in (jB6|) 

should be written as half the sum of two terms, one is as shown corresponds to and the other corresponds to 
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